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In paper \T\ with using of the Cuntz algebra representation some methods of 
construction of wavelets with scaling coefficient N >2 are considered. In paper 
[2] it is shown, a construction of wavelets at the prescribed scaling function 
(p{x). In this paper a simple method of construction of scaling function (p{x) 
and orthogonal wavelets with the compact support for any natural coefficient 
of scaling > 2 is given. Examples of construction of wavelets for coefficients 
of scaling N = 2 and = 3 are produced. 

1. Scaling functions and wavelets. Let > 2 is an integer, Z is set of all integers 
and L^(M) is Hilbert space of square integrable functions. 

Definition 1. Function <f{x) G L^(]R) is called N-scaling, if it can he represented as 

ip{x) = VNj2hn^iNx-n), (1) 

neZ 

where coefficients hn, n G Z satisfy to condition Yin l^nP < relationship (1) is 

called the N -scale equation (refinement equation). The set {hn} of coefficients of expansion 
in the equation (1) is called the scaling filter. 

Note 1. If A^-scaling function (p{x) has the compact support of length L, then the sum 
in equation (1) is finite, contained at most L[N — 1) + 1 components. 

The Fourier transform of A^-scale equation is 

f(.)^„,Q^Q, (2) 

where 

Ho{^) = A^.Y.^ne-'''^ . (3) 

The function Hq{u) is called frequency function of scaling function (p{x). 

In the orthogonal case translations of scaling function (p{x — n), n G Z form orthonormal 
basis of the subspace Vq in L^(M), and translations = \fN^p{Nx — n), n G Z on 1/iV, 

form orthonormal basis of the subspace Vi in L^(]R). Thus Vq <zVi. In the orthogonal case 
to the scaling function ^p{x) corresponds A^— 1 wavelets-functions ip^{x)..., ip'^~^{x), for each 
of which translations ipQ^{x) = ip^{x—n), n G Z form orthonormal basis of subspaces Wq in 

L^(]R), and expansion in the direct sum of orthogonal subspaces Vi = Vo®Wq(B- ■ -(BW^'^ 
be valid. 



Wavelets ip^lx)..., '^{x) form orthonormal basis L^(R): 

{'iPlSx) = ^xlj\N^x - n), J, n G Z, A; = 1, 2, . . . , iV - 1}. 

As wavelets ip^{x)..., ip^'^^x) belong to space Vi they are decomposed on basis of this 
space, 

^P'ix) = VNJ29n^iNx-n). (4) 
The coefficients {g^} is called filters of wavelets ip^{x), k = 1,2, . . . , N — 1. Let 

- the frequency functions corresponding to wavelets ip^{x)..., ip^'^i^x). The Fourier trans- 
form of equalities (4) is 



For the frequency functions Hk{uj) the following matrix is unitary [2], [2], 

/ Ho{z) Hoipz) 
^ _ Hi{z) Hi{pz) 



Hoip'^-'z) \ 
H.ip^'-'z) 



\H^-i{z) Hm-i{pz) ... Hm-i{p''~^z) J 



(6) 



where z = e '^^ and p = e ^"^^/^ , The matrix (6) has special view. It is possible to avoid 
of this special view of the matrix H{z) with Fourier transform on cyclic group Z/NZ = 
{1, p, p^ . . . , p^-^} [2]. We shall define 



(7) 



It is easy to see, that the sum on the right depends from w = z^ . Also transformation 
(7) accurate within coefficient \fN is sample of elements with degrees z^^ in polynomials 
Hk{z), z~^Hk{z)..., z~^^^Hk{z). Inverse transformation is defined by the formula [2] 

N-l 



j=0 



From last relation we shall obtained the following matrix equality: 



/ 1 



Hiz) 



A{z 



N\ 



1 

pz 



p^~'z 



\ 



A{z^)R{z). (9) 



In this expression the matrix ^4(2;^) is already arbitrary unitary matrix with polynomial 
elements. Now specificity of the matrix H{z) go to the matrix 



/ 1 



R{z) 



1 

pz 



\ 



p^'^z 



(10) 



Let's mark, that the matrix R{z) is unitary on the unit circle z = e"*"^. 

Specifying the polyphase matrix A{w), we can construct the matrix of frequency func- 
tions H{z) by the formula (8) and, together with it, frequency functions of wavelets 
Hi{z),. . . ,Hiy_i{z), hence, and wavelets ip^{x) . . . ,ip^~^{x). 

In work [2] the scheme of construction of the polyphase matrix A{z^) is given in the 
supposition, that polynomial frequency function Hq{z) is prescribed. Then it is possible to 
consider, that the first row of the matrix AQj{z^) is known. 



and it is necessary to construct remaining row of the matrix A{w). 

In the given work we shall give the simple scheme of construction of the unitary ma- 
trix A{w) which elements are polynomials with real coefficients. It allows to define both 
the scaling function (p{x) with compact support and with scaling coefficient N > 2, and 
orthogonal wavelets tp'^{x) . . . ,iIj^~'^{x). 

2. Scheme of wavelets construction. From above constructions and methods of 
work [2j follows that orthogonal systems of wavelets can be determine by the unitary matrix 
A{w) with polynomial elements with using of the formula H{z) = A{z^)R{z), where R{z) 
- the special matrix (10). We shall give the simple method of construction enough big 
set of unitary matrixes A{w) with polynomial elements. It will allow to obtain both the 
A^-scaling function with the compact support, and orthogonal wavelets. 

Let's choose any orthogonal matrix Aq = {ajj, i, j = 0, 1, . . . , A^ — 1} of the order N >2. 
We shall multiply it on the diagonal unitary matrix Dkiw) = dia.g{w^° ,w''^ , . . . ,w^^-^), 
where k = {ko, ki, . . . , kN-i) is set of integers and \w\ = 1, and then - on the orthogonal 
matrix Bq = {bij, i, j = 0,1, . . . , N — 1}. In outcome we shall obtain unitary matrix 



which elements, Aij = Yls=o ^isbsjw''% are polynomials on the variable w with real coeffi- 
cients. 

Now we shall substitute w = z^ , where z = e~^'^. We shall obtain the unitary matrix 
A{z^) with polynomial elements and real coefficients. We shall multiply it on the unitary 
matrix R{z). Then we shall obtain the unitary matrix H{z) of frequency polynomial func- 
tions Hq{z), Hi{z)..., H]\r^i{z) of orthogonal system of wavelets (p{x), 'ip^{x) . . . ,'ip^~^{x), 
where the first function (p{x) is scaling, and remaining - wavelets. Thus, 




(11) 



A{w) = AoDk{w)Bo, 



(12) 



/ Ho{z) Hoipz) ... Ho{p''-'z) \ 
H,{z) H,{pz) ... H^ip^'-'z) 



AoDk{z'')BoR{z). 



(13) 



\Hr,.i{z) H^-iipz) ... Hr^-iip^'-'z) J 
From (13) follows the expression for frequency functions: 



N -I. 



(14) 



s,j=0 



In order to the obtained the functions Hk{z) would be frequency functions of orthogonal 
wavelets, it is necessary, that the sum of coefficients for Hq{z) would be equal to unit, and 
the sums of coefficients for remaining functions Hi{z)..., Hn-i{z) would be equal to zero: 



^ iV-1 ^ N-1 

-r= Aksbsj = -7= (Iks = 0, A; = 0, 1, . . . , AT - 1. 

V-'' „ n y ly r, 



,i=o ^ .s=o j=o 

These equalities can be represented in tlic matrix view: 

/ floo C^Ol ■ ■ ■ Oo,Ar_i \ / 

aio Oil ... ai,iv-i 



aAr-i„Ar-i / 



&00 + ■ ■ ■ + ^0,Af-l 

6io H h biN-i 



\ 



\ ^iV-1,0 + • • • + &Ar-l„Ar-l / 





V I 



(15) 



Choosing various orthogonal matrixes Aq and 5o, which satisfy the equality (15), we 
obtain various frequency functions of wavelets (14). 

For construction enough simple class of orthogonal wavelets with the compact support 
and scaling coefficient N > 2, we shall take as an orthogonal matrix Aq the following 
matrix: 



1/VN 
1/V2 

1/Vq 



1/VN 
-1/V2 

1/Vq 



1/VN 


-2/V6 



1/VN 





\ l/-s/N{N -1) l/-s/N{N -1) l/-s/N{N-l) ... -{N - l)/-s/N{N - 1) J 



The matrix Aq transform vector of units e = (1, 1, . . . , 1) to the vector ^/Ncq = 
0,...,0), Aqc = VNcq. Then from equality (15) follows, that elements of the 
orthogonal matrix Bq should satisfy to the following system of equations: 



^00 + ^01 + • • • + bo^N-1 

bio + bii + ... + 6i,iv-i 



(16) 



bN-1,0 + bN-1,1 + . . . + bN-l,N-l 



The solution of this system will be any set of orthonormal vectors (rows) which coordi- 
nates satisfy to the equation of the plane Xo + Xi + . . . + X]\i^i = 1 in M . It is obvious, that 
coordinates of basis vectors Cq = (1, 0, . . . , 0), Ci = (0, 1, 0, ... , 0)...,e7v-i = (0, . . . , 0, 1) 
satisfy to this equation. The given solution corresponds to the identity matrix Bq. Any 
other solution can be obtained by rotation of the basis solution cq, ei, . . . , e^v-i around of 
vector e = Co + ei + ■ ■ ■ + e^^i, i.e. in the plane Xq + Xi + . . . + xjv_i = 1. We shall find 
these solutions. We shall take rotation around of axis Oxq: 



M 



( 1 
















1 



AT-l 



m 



N-l 



(17) 



As Aqc = VNgq, then rotation around of axis e is given by the matrix Mg = Aq^MAq. 
Then rows of the matrix Bq will consist of coordinates of vectors-columns which are obtained 
from Co, ei, . . . , ejv-i by action on them matrix Mg. Therefore the matrix Bq is transposed 
to Me. Then 

Hm(z) = AoDk(z'')MjR(z) = AoDk(z'')A'^M^AoR{z), (18) 
where M - any orthogonal matrix of view (17) and -Dfe(u') = diag(w'^", w'^'^, . . . , w''^-'-). 
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sin t cos t / ' \ — sin u cos u 



The formula (18) gives the direct method of construction the big family of frequency 
functions Hq{z), Hi{z)..., Hj^_i{z) and orthogonal wavelets with the compact support (p{x), 
ip^{x) . . . , ip^~^{x). Wavelets of the family depend of the orthogonal matrix M of view (17) 
and of the vector of degrees k = {ko, ki, , kjy^i) which it is possible to set arbitrarily. 

3. Construction of orthogonal wavelets with compact support for N = 2. 
In the given section we shall show by the example of scale N = 2 effectiveness of the 
wavelets construction scheme explained above. Though the matrix Dk{w) can be anyone, 
we shall take for example the diagonal matrix Di{w) = diag{l,w}, \w\ = 1. In case N = 2 
orthogonal matrixes Aq and Bq can be in the general view: 

■ / cost sint \ / cosM sinw 

Then 

^, \ _ ^ f cost sint \ / 1 \ / cosu sinw \ / 1 1 

y/2 \ — sin t cos t J y z"^ J y — sinu cos u J y z pz 

Frequency functions are 

Hq{z) = (cost cos u + (cost sin m)^; — (sin t sin m)^;^ + (sin t cos m) 2;^) , (19) 
v2 

Hi{z) = —z= (— sintcosM — (sintsinM)^; — (costsinu)z^ + (cost cos m)z'^) , (20) 
v2 

The sum of coefficients of frequency function Hq{z) should be equal to unit, and the 
sum of coefficients of frequency function Hi{z) should be equal to zero. The system (15) 
becomes: 

cos t sin t \ / cos -u + sin -u \ _ / -\/2 
— sin t cos t J y cos m — sin -u y \^ 

cos M + sin M = a/2 cos t 
cos u — sin M = \/2 sin t 

Solving last system, we obtain, u = tt/A — t. 

Thus, we have constructed the family of frequency functions of the wavelets specified 
by formulas (19), (20) in which u = n/A — t. After elimination of the variable u, we obtain:: 

Ho{z) = i (1 + cos2t + sin2t + (1 + cos2t - sin2t)2;+ 
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+ (1 - cos2t - sin2t)z^ + (1 - cos2t + sin2t)z^) , (21) 



Hi{z) = i (-1 + cos2t - sin2t + (1 - cos2t - sin2t)z+ 

+ (-1 - cos 2t + sin 2t)z^ + (1 + cos 2t + sin 2t);z^) . (22) 

The given frequency functions Hq{z) and Hi{z) coincide with the same, but obtained 
other methods in work [T] . Various wavelets of Haar, Daubechies wavelets and their analogs 
include into this family. In the following section some examples are given. 

Choosing other matrix Dk{z'^), similarly we can construct other orthogonal wavelets 
with other support length. 

4. Examples of scaling functions and wavelets for N = 2. We shall calculate 
values of coefficients of the obtained frequency functions (19), (20) for various parameters 
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t and u and we shall find corresponding filters and wavelets (fi{x) and ip{x). From formulas 
(21), (22) follows what enough to take parameter values t on interval of length tt. We shall 
consider the following parameter values t: 0, ±7r/12, ±7r/6, ±7r/4, ±7r/3, ±57r/12, tt/2. 
4.1 Parameter values i = 0, u — n/A. Coefficients of wavelets filters: 

ho = -^(1, 1, 0, 0), 51 = -^(0, 0, -1, 1). 

We have obtained wavelets of Haar with the support on unit interval. Refinement equations: 
ip{x) = ip{2x) + ip{2x - 1) and if^ix) = -^p{2x - 2) + ^l){2x - 3). 

4.2. Pcirameter values t — 7r/4, u — Q. Coefficients of wavelets filters: 

ho = -^(1, 0, 0, 1), 51 = -^(-1, 0, 0, 1). 

We have obtained wavelets of Haar with the support on interval [0,3]. Refinement equa- 
tions: ip{x) = <f{2x) + ip{2x — 3) and ijj{x) — — -0(2^) + -0(20; — 3). 

4.3. Pcirameter values t — 7r/2, u — 7r/4. Coefficients of wavelets filters: 

ho = -^(0, 0, 1, 1), 51 = -^(-1, 1, 0, 0). 

This is wavelets of Haar. Scaling function has the support on interval [2,3]. Refinement 
equations: (p{x) = ip{2x — 2) + ip{2x — 3) and ip{x) = —ip{2x) + tp{2x — 1). 

4.4. Pcirameter values t — 7r/4, u — 7r/2. Coefficients of wavelets filters: 

ho = -^(0, 1, 1, 0), 51 = -^(0, 1, -1, 0). 

This is wavelets of Haar. Scaling function has the support on interval [1,2]. Refinement 
equations: ip{x) = ip{2x — 1) + (p{2x — 2) and iIj{x) = 'il){2x — 1) — '0(2a; — 2). 

4.5. Pcirameter values t — 7r/12, u — 7r/6. Coefficients of wavelets filters: 

^0 = ^(3 + V^,l + V^, 1-^,3- V^), 51 = ^(-3 + V^,l-^,-l-^,3 + V^). 

o o 

The result will be wavelets with coefficients which are obtained by permutation of co- 
efficients of the classical Daubechies wavelets with the support of length 3. Refinement 
equations: 

3 + ^3,, 1 + ^3,^ , 1-V^, , 3-^3, 
ip{x) = -^^{2x) + -^v{2x - 1) + -^^{2x - 2) + -^v{2x - 3). 

In figure 1 graphs of wavelets are shown. 

4.6. Pcirameter values t — 57r/12, u — — 7r/6. Coefficients of wavelets filters: 

/2 /2 
/io = ^(3-^3,1 -v^,l + v^, 3 + ^3), 51 = ^(-3-v^,l + V3,-l + v^,3-v^). 

o o 

This example differs from previous only that coefficients of the filter go upside-down. 
In this case scaling function can be obtained from scaling function of example 4.5 with the 
using of argument replacement: ip{'i — x). It follows from the fact: if y?(x) - scaling function 
with the compact support [0, L] and the filter then function ip{L — x) also is scaling 
with the filter {/ii_„}.The corresponding wavelet also can be obtained from previous as: 
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Figure 1: Graphs of functions (p{x) and ip{x) for t = 7r/12, u = tt/Q 





Figure 2: Graphs of functions ip{x) and ■0(2^) for t — 57r/12, = — 7r/6 



— ■0(3 — x). The graph of scahng function (p{x) can be obtained from the graph of the Fig.l 
by mirroring about the hne x = 3/2. For the graph of the wavelet ip{x) it is necessary to 
add still mirroring about axis Ox (Fig. 2). 

4.7. Pcirameter values t — — 7r/12, u — n/S. Coefficients of wavelets filters: 



V2, 



V2, 



ho = ^{1 + V3,3 + V3,3-V3,1-V3), gi = -V(-l + v^,3-v^, -3-^3, 1 + ^3). 

o o 
The result will be Daubechies wavelets with the support of length 3. Refinement equation: 

l + \/3 , , 3 + v^ , , 3-^ , , 1-^ , 
ifix) = -^^{2x) + -^^{2x - 1) + -^^{2x - 2) + -^ip{2x - 3). 

In figure 3 graphs of wavelets arc shown. 

4.8. Pcirameter values t — — 57r/12, u — 27r/3. Coefficients of wavelets ffiters: 

/2 /2 
/io = ^(l-\/3,3-\/3,3+\/3,l + \/3), gi = ^(-1 - ^3, 3 + Vs, -3 + \/3, 1 - \/3). 

o o 

This example differs from the previous only that coefficients of the filter {hn} go upside- 
down. In this case scaling function can be obtained from Daubechies scaling function with 
the help of argument replacement: ip{3 — x), and wavelet is —ip{3 — x). 

4.9. Parameter values t = tt/G, u = 7r/12. Coefficients of wavelets ffiters: 



ho = ^(3 + \/3,3-\/3,l-\/3,l + \/3), 
o 



gi 



V2^ 



-1 - \/3, 1 - \/3, -3 + v^, 3 + V3) . 





Figure 3: Graphs of functions ip{x) and ip{x) for t = — 7r/12, u = tt/S 
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Figure 4: Graphs of functions (p{x) and 'ip{x) for t — n/Q, u — t:/12 

The result will be wavelets with coefficients which are obtained by permutation of Daubechies 
wavelets coefficients. The refinement equation: 

if{x) = — ^ — ip{2x) + — ^ — ip{2x - 1) + — ^ — ^p{2x - 2) + — ^ — ip[2x - 3). 

In figure 4 graphs of wavelets are shown. 

4.10. Parameter values t — tt/S, u — — 7r/12. Coefficients of wavelets filters: 

/2 /2 
/io = ^(l + \/3,l-\/3,3-\/3,3 + \/3), ^1 = ^(-3-V3,3-\/3,-l + \/3,l + \/3). 

o o 

This example differs from the previous only that coefficients of the filter go upside- 
down. In this case scaling function can be obtained from the previous scaling function by 
replacement of argument: (^(3 — x), and wavelet is — ^^(3 — x). 

4.11. Parameter values t — — 7r/3, u = 77r/12. Coefficients of wavelets filters: 

/2 /2 
/io = ^(1-a/3,1 + \^,3 + V3,3-a/3), = ^(-3 + a/3,3 + a/3, -1-a/3, 1-^). 

o o 

The result will be wavelets with coefficients which are obtained by coefficients permutation 
of Daubechies wavelets with the support of length 3. Refinement equations: 

^[x) = -^(p{2x) + -r^(p{2x - 1) + — -r^(p{2x - 2) + -r^^{2x - 3). 



Figure 5: Graphs of functions (f{x) and ip{x) for t — 



-7r/3, u = 77r/12 



In figure 5 graphs of wavelets are shown. 

4.12. Parameter values t — — 7r/6, u = 57r/12. Coefficients of wavelets filters: 

ho = ^(3-\/3,3 + \/3,l + \/3,l-\/3), gi = ^(-1 + \/3, 1 + \/3, -3- Vs, 3- v^). 
8 8 

This example differs from the previous only that coefficients of the filter go upside- 
down. In this case scaling function and wavelet can be obtained from the previous by 
replacement of argument: ip{3 — x), —ip{3 — x). 

5. Construction of wavelets in case N — 3. In this section we shall show the scheme 
of scaling function and wavelets construction for = 3. Though the diagonal matrix 
Dkiw) can be anyone, we shall take for example the diagonal matrix Di{w) = diag(l, w, 1), 
\w\ = 1. The matrix Aq is: 

1/^3 l/VS l/VS 

1/V2 -1/V2 

l/y/6 1/V6 -2/\/6 
Elements of the second orthogonal matrix Bq should satisfy to conditions: 
boo + boi + bo2 = 1, bio + bu + 612 = 1, &20 + &21 + ^22 = 1- 

The solution of this system will be any set of orthonormal vectors which coordinates satisfy 
to the equation of the plane -I-X2 = 1. It is obvious, that coordinates of basis vectors 

ei, 62, 63 satisfy to this equation of plane. For this solution the matrix Bq it is identity. 
And we obtain the wavelets of Haar, 

/ 1/Vs w/Vs 1/Vs \ 

AoDi{w)Bo = I/V2 -w/V2 , (23) 
V 1/V6 w/^/^ -2/^ / 

/ 1/V3 zVVs l/x/3 \ / 1 1 1 
H{z) = -= -zyV2 \\zpzp'z 

1/a/6 zyVe -2/V6 J \ p'z' p'z^ 

Ho{z) = ^(1 + z' + z% Hi{z) = i=(l - z'), H,{z) = -i=(l - 2^2 + z% 

The maximum degree of frequency function Ho{z) is equal to four, the support length 
L is equal to two, as it is find from the formula L{N — 1) + 1 = deg{Ho{z)) + 1. 
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It is easy to see, that scaling function if{x) is characteristic function of interval [0,2), 
(p{x) = Xlo,2){x). The refinement equation and wavelets (Fig. [6]): 

(p{x) = (p{3x) + Lp{3x — 2) + {p{3x — 4), 
^^(x) = ^(^(3a:)-^(3x-4)), 

^P'ix) = (ifiSx) - 2ipi3x - 2) + ifiSx - 4)) , 

Any other solution can be obtained by rotation of basis vectors Cq, ei, 62 in the plane 
xo + 2;i + 0:2 = 1. We shall find these solutions. Let 

10 
M(t) = I cost -sint 
sin t cos t 

- the matrix of rotations around of the axis Cq. Then 

T / 1 + 2 cos t 1 — cos sin t 1 — cos if: — sin t 

Me{t) = A^^M{t)Ao = - \ 1 - cost - v^sint l + 2cost 1 - cost + ^3 sint 



3 



1 — cos t + -\/3 sin t 1 — cost — sin t l + 2cost 



Let's make rotation Me(t)ek of column vectors cq = (1, 0, 0), ei = (0, 1, 0), 62 = (0, 0, 1), 
and we obtain rows of the required matrix BQ{t): 

^ / 1 + 2 cos t 1 — cos t — y/S sin t 1 — cos t + \/3 sin t 
^o(^) = o I 1 -cost + VSsint l + 2cost 1 - cost - \/3sint 1. (24) 
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1 — cos t — -\/3 sin t 1 — cost + sin t l + 2cost 



Then H(t,w) = AQDi{w)Bo{t)R{z) where the matrix AqDi^z'^) is represented by the 
formula (23), -Bo(''^) ^ by the formula (24) and the matrix R{z) is 

^,11 1 

R{z) = —= z pz p^z 



z^ p^z^ p'^z^ 



Multiplying all these matrixes and choosing elements of the first column, we obtain, 

Ho{t,z) = - (2 + cost - v^sint + (2 - 2cost)z+ (2 + cost + v^sint)2;^+ 
9 V 

+ (1 - cost + sin t)z^ + (1 + 2 cost)/ + (1 - cost - sin t)/j , (25) 

Hi(t, z) = — ^ (1 + 2 cost + (1 — cost — y/3smt)z + (1 — cost + sint)z^ — 
3v6 ^ 

-(1 - cost + VSsint)/ - (1 + 2 cost)/ + (-1 + cost + V3smt)z^^ , (26) 
H2{t,z) = (^-1 +4cost + 2v^sint- (1 - cost + 3^3 sin t)z- 
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Figure 6: Graphs of functions ip{x), ip-^lx) and ip'^{x) for t = 

-(1 + 5 cost + A/3sint)z^+ 
+ (1 - cost + V3smt)z^ + (1 + 2 cost)/ + (1 - cost - VSsint)/ 



(27) 



6. Examples of scaling functions and wavelets for = 3. We shall calculate 
coefficients of the obtained frequency functions (25), (26) and (27) for various parameter 
values t. The obtained filters allow to find corresponding wavelets v?(x), ip^{x) and ip^{x) 
by usual methods [5], [3]. It is enough to find scaling function f{x). Wavelets - functions 
ip^{x) and ip^{x) are defined by formulas 

iIj\x) = VNj^dn^iNx-n), il;\x) = VNj^dMNx-n) 



with known filters {Qn} and {g^} and function ip{x). 

Let's consider the following parameter values t: 0, n/6, n/A, vr/S, 7r/2, 27r/3, vr, 47r/3. 
For each case graphs of wavelets-functions are shown. 

6.1. Value of parameter t = 0. This case has already been considered above. It is 
wavelets of Haar with the support [0, 2] (Fig. [6]). 

6.2. Value of parameter t = tt/Q. Filters of scaling function (p{x) and wavelets ip^lx) 
and ip'^{x): 

/3 

/lo = ^(2, 2 - VS, 2 + VS, 1, 1 + ^3, 1 - ^3), 
9 



9i 



V6 



(3 + VS, -3 + VS, VS, - V3, -3 - V3, 3 - VS), 



^2 



76 



-1 + 3^3, -1 - VS, -1 - 2^3, 1, 1 + VS, 1 - VS). 



The refinement equation: 

^(x) = \{2^{3x) + (2 - V3)<^{3x - 1) + (2 + V3)^{3x - 2) + (^(32; - 3) + 

+ (1 + V3)^i3x - 4) + (1 - V3)ip{3x - 5). 
Graphs of wavelets are shown in figure [7] 

6.3. Value of parameter t = 7r/4. Filters of scaling function (p{x) and wavelets iIj^{x) 
and ip'^{x): 

\/3 

/lo = (4 + 72 - Ve, 4 - 2^2, 4 + 72 + 76,2-^2 + ^6,2 + 2^2, 2-72-76), 

io 
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Figure 7: Graphs of functions (f{x), ip-^lx) and il''^{x) for t = 7i/6 



Figure 8 
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: Graphs of functions (f{x), ip^lx) and ip'^{x) for t = 7r/4 



c/i = — (2^3 + 2^6, -3^2 + 2^3 - ^6, 3^2 + 2^3 - ^6, -3^2 - 2^3 + ^6, 
36 



-2^3 - 2^6, 3^2 - 2^3 + ^6), 



\/6, 



^2 = — i-2 + + 2^6, -2 + - 3v^, -2 - 5v^ + v^, 2 - + v^, 
3o 

2 + 2^2,2-72-^6). 

Graphs of wavelets are shown in figure [8] 

6.4. Value of parameter t = 7r/3. Filters of scaling function ^{x) and wavelets ^^(x) 
and il)'^{x) are 

/io = ^(1,1,4,2,2,-1), 
= ^(2, -1, 2, -2, -2, 1), g2 = ^(4, -5, -2, 2, 2, -1). 



6 

The refinement equation: 



1^ 



(^(x) = -(<^(3x) + (^(32; - 1) + 4(^(32; - 2) + 2^(32; - 3) + 2(^(3x - 4) - ^{"ix - 5). 
Graphs of wavelets are shown in figure [H 

6.5. Value of parameter t = tt/2. Filters of scaling function ^p{x) and wavelets il)^{x) 
and ip'^i^x): 



12 





Figure 9: Graphs of functions (f{x), ip-^lx) and ip'^lx) for t = 7i/3 
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Figure 10: Graphs of functions (f{x), ip-^lx) and il''^{x) for t = 7r/2 



/lo = ^(2 - 73, 2, 2 + V3, 1 + v^, 1, 1 - 73), 
9 

/fi 

91 = T^(73, -3 + 73, 3 + 73, -3 - 73, -73, 3-73), 
18 

76 

92 = + 273, -1 - 373, -1 + 73, 1 + 73, 1, 1 - 73). 

lo 

The refinement equation: 

V9(x) = ^((2 - VsMSx) + 2y;(3x - 1) + (2 + 73)^(3x - 2) + (1 + Vs^Sx - 3) + 

+^{3x - 4) + (1 - V3)ip{3x - 5). 
Graphs of wavelets are shown in figure dUl 

6.6. Value of parameter t = 27r/3. Filters of scaling function (p{x) and wavelets 
ip^i^x) and ip'^i^x) is 

/io = -^(0,1, 1,1, 0,0), 
g, = -^(0, 0, 1, -1, 0, 0), 92 = -^(0, -2, 1, 1, 0, 0). 
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Figure 11: Graphs of functions (f{x), ip-^lx) and il''^{x) for t = 27r/3 
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Figure 12: Graphs of functions (f{x), ip^{x) and ip'^{x) for t = n 



It is wavelets of Haar. The scahng function (p{x) is characteristic function of interval 
[1/2,3/2), (p{x) = X[i/2,3/2)(a^)- The refinement equation and wavelets: 

(p{x) = ip{3x — 1) + ip{3x — 2) + {p{3x — 3), 

^i(x) = ^(^(3x-2)-v^(3x-3)), 

tp^ix) = ^(-2^(3x - 1) + y?(3x - 2) + ip{3x - 3). 
v2 

In figure [TT] graphs of wavelets are shown. 

6.7. Value of parameter t = tt. Filters of scaling function (f{x) and wavelets ip^lx) 
and ip'^{x): 

\/3 

/io = ^(1,4,1,2,-1,2), 

91 = ^(-1, 2, 2, -2, 1, -2), g2 = ^(-5, -2, 4, 2, -1,2). 
o 18 

The refinement equation: 

^(x) = \{^{3x) + 4:<^{3x - 1) + ^{3x - 2) + 2^{3x - 3) - (^(32; - 4) + 2^{3x - 5). 



Graphs of wavelets are shown in figure 
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Figure 13: Graph of function (f{x) for t = Att/3 



6.8. Value of parameter t = An/S. Filters of scaling function v^(x) and wavelets 
ip^{x) and ip'^{x): 

/zo = -^(1,1, 0,0, 0,1), 

g, = i=(0, 1, 0, 0, 0, -1), 92 = -^(-2, 1, 0, 0, 0, 1). 
The refinement equation and frequency functions: 

(p{x) = ip{3x) + ip{3x — 1) + f{3x — 5), 

Ho{z) = ^(1 + ^ + ^'), ^i(^) = - ^')' ^2 W = ^(-2 + z + z'). 

Let's mark, that scaling function (f{x) has a complicated structure. Its support has fractal 
properties. 

References 

[1] O. Bratteli, D. E. Evans, and P. E. T. Jorgensen, Compactly supported wavelets and 
representations of the Cuntz relations. Appl. Comput. Harmon. Anal. Vol. 8, 2000, 
166-196. arXiv.org:math.FA/9912129 

[2] O. Bratelli, P. E. T. Jorgensen, Wavelet filters and infinite-dimensional unitary groups. 
arXiv.org:math.FA/0001171 v3, (2000), 31 p. 

[3] I. Daubechies, Ten Lectures on Wavelets. CBMS-NSF Regional Conf. Ser. in Appl. 
Math., vol. 61, Society for Industrial and Applied Mathematics, Philadelphia, 1992. 

[4] P. N. Podkur, N. K. Smolentsev, Construction of some types wavelets with coefficient 
of scaling N. arXiv.org:math.FA/0612573, 2006, 19 P. 

[5] N. K. Smolentsev, Osnovy teorii weivletov. Weivlety v MATLAB. M., DMK Press, 2005 
(in Russian). 



15 



